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Abstract 



I review the status of the modern theoretical approach to weak decays of heavy 
flavor hadrons based on the 1/tuq expansion in QCD. The qualitative features are ex- 
plained and the subtleties in simultaneously incorporating perturbative and power- 
suppressed effects are addressed. A few topical phenomenological applications are 
discussed in quantitative detail. 



"Invited talk at the Symposium on Radiative Corrections (CRAD 96), Cracow, 1-5 August 1996. 



1 Introduction 



The Standard Model (SM) of fundamental interactions has so far yielded a consistent 
description of various experimental data; still, from the general perspective it is often 
viewed incomplete. A key role in exploring the SM is played by studying electroweak 
decays of heavy flavor particles. Unique information about the masses and mixing 
of quarks is obtained here thus setting the stage for addressing some of the most 
mysterious questions about the origin of the quark mixing and CP violation. 

The fact that the dynamics of heavy flavor hadrons represents a peculiar field 
where the strong interaction effects allow a systematic treatment based on first 
principles of QCD, was realized already in the 70s. Yet, the full variety of theoretical 
tools available for gauge theories was applied here only over the last few years. 

Two main directions were instrumental in developing QCD towards its current 
status of the theory underlying strong interactions: one was based on studying sym- 
metry properties like isotopic or SU (3) symmetry and chiral invariance; the second 
exploited the asymptotic freedom which allowed using dynamical calculations in the 
perturbative way for high-energy processes. These two approaches, the 'symmetry- 
based' and 'dynamical', are clearly seen in the evolution of the heavy quark theory. 

During the early period to the end of the 80's, most applications carried the 
'dynamical' spirit and were done often at a somewhat simplified intuitive level. The 
nonperturbative effects typically were thought to be small even in the decays of 
charm particles. At the same time, the basic facts about the heavy quark spin and 
flavor symmetry were realized and incorporated to the necessary extent. 

The symmetry considerations flourished for a few years with the formulation of 
the so-called Heavy Quark Effective Theory (HQET) which made these symmetries 
explicit at the level of the Lagrangian. It provided a convenient framework for dis- 
cussing exclusive semileptonic or electroweak radiative decays and calculating basic 
perturbative corrections to them. On the other hand, the range of applications of 
HQET was limited, and not all virtues of the general Heavy Quark QCD Expansion 
(HQE) were employed. 

Finally, over the last few years a consistent well-defined dynamical approach 
has been developed, which automatically respects the heavy quark symmetries in a 
manifest way. It has been realized that the nonperturbative corrections are often 
sizable even in the decays of beauty. Incorporating the actual dynamical calculations 
allowed one to make the most precise determinations of \V c b\ and \ V u b\. 

The main effects in a weak decay of a heavy quark Q originate from distances 
~ 1/mq <C 1/Aqcd- Belonging to the weak coupling domain, a s (mg) <C 1, they 
are tractable through perturbation theory. The QCD-interaction becomes strong 
only when the momentum transfer is much smaller than the heavy quark mass, 
k <C itlq. This simple observation elucidates the two basic ingredients of the heavy 
quark expansion that is genuinely based on QCD: 

• The nonrelativistic expansion, which yields the effects of 'soft' physics in the 
form of a power series in 1/toq. 
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• The treatment of the strong interaction domain based on the Operator Product 
Expansion (OPE). 

Unless an analytic solution of QCD is at hand, these two elements appear to be 
indispensable for the heavy quark theory. 

The necessity of separating the two domains characterized by momenta much 
larger than Aqcd and those which do not scale with rriQ, is thus manifest. The 
perturbative corrections, on the other hand, bridge the short-distance domain with 
the least understood domain of confining strong forces. 

The general idea of separating two domains and applying different theoretical 
tools to them was formulated long ago by K. Wilson [0J in the context of problems in 
statistical mechanics; in the modern language, applied to QCD it is similar to lattice 
gauge theories. The treatment of essentially Minkowski quantities has, however, 
some peculiarities which attracted attention of theorists very recently. 

The perturbative corrections - in spite of being often complicated technically - 
look straightforward from the conceptual viewpoint. On the contrary, getting even 
a minimal model-independent information about the nonperturbative physics often 
looks like a miracle. In view of the nature of this Symposium I will still try to 
put more accent onto the perturbative aspect, in particular, since the perturbative 
corrections got a renewed attention during the last two years. The overlap between 
the 'perturbative' and 'nonperturbative' effects is non-trivial indeed. 

There are different technical ways to perform I/tuq expansion in QCD; the most 
popular schemes are known as the Nonrelativistic QCD (NRQCD) || ||| and the 
Heavy Quark Effective Theory (HQET) . In what concerns the treatment of the 
perturbative corrections, in the standard HQET they are simply added to what is 
considered as the nonperturbative physics. On the other hand, the application of 
Wilson's idea requires separating of any observable quantity into short- distance and 
long-distance pieces which, formally, is a different procedure. Schematically, the 
HQET decomposition 

An observable = perturbative piece + nonperturbative piece 

is replaced in Wilson's OPE by a separation-scale-dependent representation 

An observable = short-distance piece + long-distance piece . 

The novel feature we face in the theoretical analysis of beauty decays nowadays 
is that they often allow - and even demand by virtue of the existence of precise 
experimental measurements - rather accurate predictions, requiring a simultane- 
ous treatment of perturbative and nonperturbative QCD effects with the enough 
precision in both. This problem is not new; the theoretical framework has been 
elaborated more than 10 years ago ||, but its phenomenological implementation 
was not mandatory until recently. It is the exceptional beauty of the fifth quark 
that allows - in spite of interference of the nonperturbative effects - to make cer- 
tain reliable theoretical predictions with errors at a percent level, which necessitates 
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taking care of such 'puristic' theoretical subtleties. Nevertheless, failure to incorpo- 
rate them properly, leads to certain theoretical paradoxes and, unfortunately, some 
superfluous controversy in the numerical estimates surfacing every now and then in 
the literature. 

In spite of the significant progress made over the last years, the theory of the 
heavy flavors is not yet a completed field and is undergoing to further extensive 
development. There are quite a few topics which separately deserve a proper review: 

• Sum rules for the weak transitions with the nonperturbative effects included. 
In a certain way they replace the quantum mechanical (QM) description in 
the field theory. The perturbative corrections to the sum rules were actively 
discussed, often with quite different conclusions. 

• The so-called BLM-resummation and renormalons. 

• The QCD-based OPE vs. 'Fermi motion' in heavy quark decays. 

• The QCD sum rules calculations of the most interesting electroweak formfac- 
tors in the b — > u (b — > s) transitions. 

• Attempts to obtain constraints on the transition formfactors relying on the 
most general unitarity and analyticity properties. Without an actual dy- 
namical input they are too weak, however. There is an interesting idea to 
incorporate additional constraints from lattice QCD. 

• The heavy quark expansion for the production of heavy flavor hadrons and 
quarkonia. 

The practical applications were in focus, too: 

• Determination of jV^,] and \V u b\. 

• A model-independent extraction of m h . 

• Radiative corrections to BR sl (_B). 

There is also a renewed interest in the inclusive nonleptonic widths. 

Finally, the violations of duality is more and more discussed in the context of 
the heavy quark physics. 

Even this, by necessity truncated, list gives a feeling of the extensive develop- 
ment taking place. Inevitably, I cannot cover here all, and even a large part of 
recent theoretical advances, or simply list all relevant references. A true measure 
of our understanding of the strong interaction dynamics is, eventually, determined 
by the theoretical accuracy with which we can extract the fundamental underlying 
parameters like the KM mixing angles from experimental data. Therefore, I will 
focus in this talk on a few selected topics that illustrate the theoretical framework 
used in the most accurate determination of the flavor mixing parameters, and review 
the overall status of the heavy quark expansion for beauty particles, with the main 
emphasis on the qualitative features. 
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a b 

Figure 1: Exclusive B — > (a) and generic (b) semileptonic decays. 

2 Semileptonic decays 

The QCD-based heavy quark expansion can in principle equally be applied to all 
types of heavy flavor transitions. Semileptonic decays, however, are the simplest 
case and I shall devote most of the attention to them. For practical reasons I focus 
on the b — > c transitions; a brief discussion of the b — > u decays will be given later. 

A typical semileptonic decay is schematically shown in Figs. 1. Generally, two 
types of the decay rates can be singled out: the inclusive widths, where any com- 
bination of hadrons is allowed in the final state, and the exclusive decays, when a 
transition into a particular charmed hadron is considered, usually D or D*. 



2.1 Inclusive semileptonic widths 

The four-fermion semileptonic decay width of a heavy quark has the form 

G\m\ 2 / m\\ 



where z is the known phase space suppression factor and se generically includes all 
QCD corrections. In the heavy quark limit the difference resulting from using the 
quark mass m b and the meson mass Mb in eq. (P disappears: 

M B - m b 1 ^ 



m b m b 

For the actual b quark ml and M B differ by a factor of 1.5-2, which formally consti- 
tutes a power-suppressed nonperturbative effect. This demonstrates the necessity of 
a systematic control of the nonperturbative corrections even in the decays of beauty 
particles. 

The central result for the inclusive decay widths in QCD which is guaranteed 
by the application of the OPE to the heavy quarks in full QCD, is the absence of 
the 1/rriQ corrections |6], [7]]. This seems non-trivial a priori since there are such 
corrections to the masses, in particular, differentiating masses of different types of 
hadrons (A&, B, B**). The physical reason behind this fact is the conservation 
of the color flow in QCD, that ensures the cancellation of the effects of the color 
charge (Coulomb) interaction in the initial and final states. In terms of nonrelativis- 
tic quantum mechanics it is the cancellation between the phase-space suppression 



4 



caused by the Coulomb binding energy in the initial state, and the Coulomb dis- 
tortion of the final state quark wavefunctions. The inclusive nature of the total 
widths leads to the fact that the width is sensitive only to the strong interaction on 
the time scale ~ l/AE ~ I/mi, (AE denotes the energy release). The final state 
interaction effect is thus not determined by the actual behavior of the strong forces 
at typical hadronic distances, but only by the potential in the close vicinity of the 
heavy quark; the cancellation occurs universally, whether or not a nonrelativistic 
QM description is applicable. 

The leading power corrections start with terms 1 /m 2 ; they have been calculated 
in Ref. || [7| and are expressed in terms of the expectation values of two operators 
of dimension 5, which in the language of QM are interpreted as the strength of 
the chromomagnetic field at the position of the heavy quark and the square of the 
spacelike momentum of the heavy quark jiggling in the rest-frame of the hadron, 
respectively: 

/4 * ^{B^l^CT^B) " ( B \*» ■ 9sH g (0)\B) ~ ~(M 2 ,-M 2 ) ~ 0.35 GeV 2 

(3) 

/4 * ^W(^) 2 &|5> - (B\P 2 \B) (4) 

(the operators depend on the normalization point, but for simplicity I do not indicate 
this fact). The size of /i^. for B mesons is derived from the observed hyperfine mass 
splitting between B and B*. 

The value of fi 2 is not yet known directly; a model-independent lower bound was 
established in || [| 

/4 > /4 , ( 5 ) 
which puts an essential constraint on its possible values. This bound is in agreement 
with QCD sum rule calculations |TIJ, yielding a value about 0.5 GeV 2 and with an 
estimate [11] relying on the measured slope of the Isgur-Wise function. In the 



absence of the perturbative gluon corrections, as it happens, for example, in simple 
QM models, the expectation value of the kinetic operator /i 2 would coincide with 
the HQET parameter — Ai; they are different, however, in the actual field theory, 
where both /z 2 and [Aq depend on the normalization point. 

Including the nonperturbative corrections, the semileptonic width has the fol- 
lowing form [|, ^ |TJ, 0: 



G%ml , 2 J / /i 2 -/4\ (. m 2 c \ /4 2 a, (i) 



where ellipses stands for higher order perturbative and/or power corrections; z's 
are known phase space factors depending on m 2 c jm\. Irrespectively of the exact 
value of /i 2 , 1/mf corrections to T sl are rather small, about —5%, thus leading to 
the increase in the extracted value of \V c b\ by 2.5%; the impact of the higher order 
power corrections is negligible. 
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Good control of the QCD effects in the inclusive semileptonic widths provides 
the most accurate direct way to determine \V c b\ in a truly model- independent way. 
This method sometimes faces a traditional scepticism: which numerical value must 
be used for ra& and m c ? It appears that this practical problem has deep roots; failure 
to understand them is the major source of controversy about the heavy quark masses 
and inclusive widths often found in the literature. It will be briefly discussed below. 

In reality the precise value of is not too important, since the b —>■ c width 
depends to a large extent on the difference mj — m c rather than on mj itself; the 
former is constrained in the heavy quark expansion 

M B + 3Mb* M d + 3M D * 2 f 1 1 



ms _ mc = _ + _ , + ... M 3.50 GeV . (7) 



It also independently enters lepton spectra in semileptonic decays and can be 



directly extracted from the data [14]]. Numerically [151 , a change in m& by 50 MeV 
leads only to a 1% shift in |V^|. 

Heavy quark masses 

An accurate knowledge of the heavy quark mass becomes important if a few percent 
precision in \V c b\ is targeted. However, rather different estimates of rrib can be found 
in the literature. The reason is that the most popular scheme where certain features 
of the HQE were implemented, the HQET was based on the so-called pole mass of the 
heavy quarks. Not only was it a starting parameter of the HQET-based expansions, 
it is this pole mass that one always attempted to extract from the experimental 
data. It turns out, however, that the pole mass of the heavy quark is not a physical 
notion and cannot be even defined with the necessary accuracy in a motivated way: 



it suffers from an irreducible intrinsic theoretical uncertainty of order Aqcd [|16 
This applies therefore to the parameter A = Mb — m£ ole as well. 

At first sight this looks paradoxical and counter-intuitive: for the pole mass 
is a usual mass of a particle; for example, the value of m e quoted in the tables 
of physical constants is just the pole mass of the electron. In QCD there is no 
'free heavy quark' particle in the physical spectrum, and its pole mass is not well 
defined. The problems facing the possibilities to extract the pole mass from typical 
measurements were illustrated in Refs. |T7| and [ |T%] . 

The physical origin of the uncertainty brn^^ ~ Aqcd is the gluon Coulomb self- 
energy of the static colored particle. The energy stored in the chromoelectric field 
inside a sphere of radius R ^> 1/ tjiq is given by 

5-Ecouiomb(i?) oc / El d 3 x oc const - ^f-j- -L . (8) 

Jl/m b ~\x\<R IT R 

The pole mass assumes that all energy associated with the color source is counted, i.e. 
R — ► oo. Since in QCD the interaction becomes strong at Ro ~ 1/Aqcd, the domain 
outside Rq would yield an uncontrollable (and physically senseless) contribution to 
the mass ~ A QCD [16|] . 
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Figure 2: Perturbative diagrams leading to the IR renormalon uncertainty in wlq 6 
of the order of Aqcd- The contribution of the gluon momenta below uiq expresses 
the classical Coulomb self-energy of the colored particle. The number of bubble 
insertions into the gluon propagator can be arbitrary 



Being a classical effect originating at a momentum scale well below tuq, the 
appearance of this uncertainty can be traced in the usual perturbation theory, where 
it manifests itself in higher orders as a so-called 1/mg infrared (IR) renormalon 



singularity in the perturbative series for the pole mass [Tj| [20| . 

To model the effect of the running coupling and emergence of the nonperturbative 
low-momentum domain one can use the running coupling a s (k 2 ) in the one- loop 
diagram for the correction to the on-shell mass of the heavy quark, see Fig. 2. In 
the nonrelativistic limit k <C mq the expression simplifies and reads 

4 r d A k Ana s (k 2 ) 1 4 r d 3 k a s {k 2 ) 

6mQ ~ ~3 J J^yl^^h = si ^^r- (9) 

Expanding running a s (k 2 ) in terms of a short-distance a s (/i 2 ) one obtains for SrriQ 
the whole series in a s (fi 2 ). The coefficients grow factorially: 

5m Q 4a s (ba s \ n 11 2 

ttiq 3 7i \ 2tt J 6 6 

and the nonperturbative regime manifests itself as the uncertainty in defining a sum 
of such a series. Numerically, the calculation of the contribution of the domain 
below \i ~ 1 GeV would then look as follows: 

m P° le = m 6 (lGeV) + 5m peit ~ 

4.55 GeV + 0.25 GeV + 0.22 GeV + 0.38 GeV +1 GeV + 3.3 GeV +14 GeV + ... (11) 

Clearly, such a possibility does not look encouraging if one really needs to know the 
heavy quark mass accurately enough! 

The first two terms in eq. (0) constituting approximately 4.8 GeV can be called - 



with some reservations [17j - the 'one-loop' pole mass that enters typical calculations 



when order-a s corrections are computed. Recently, the third term was incorporated 



[21 1 in the analysis of the lepton spectra in the semileptonic 6-decays and the value 

A = M B -m p b ole ~ (390±110)MeV (12) 
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was deduced. The quoted error, according to jnj], includes "only experimental 
uncertainty" - however, such a numerical conclusion can hardly be taken sensibly. 

In spite of this irreducible uncertainty in mg le and A, the inclusive widths, 
though dependent on the mass, can be theoretically calculated since they are gov- 
erned not by the pole masses, but well-defined short-distance running masses rriQ^) 
with the Coulomb energy originating from distances > 1//J, peeled off. Moreover, 
it is precisely this short-distance running mass rather than the pole mass that can 
be extracted from experiment with, in principle, unlimited accuracy: the pole mass 
does not enter any genuine short- distance observable at the level of nonperturbative 
corrections, which can be calculated by means of the OPE [TJ|. It is worth noting 
that the irrelevance of the pole mass roots deeper than merely the problem of the 
IR renormalon; even if there were no problem in defining the pole mass in QCD, 
its infrared part distinguishing it from the short- distance mass is still foreign to the 
OPE. 

Applied to the inclusive widths, this observation suggests certain information 
about the importance of higher order perturbative corrections: if masses entering 
eq. (P) are the pole masses, the perturbative series 



rp ert = r o8e p ert = r 1 + oj^ + aa ( ^ ) + ... (13) 




is poorly behaved, with coefficients factorially growing, which makes it non- 
convergent and the radiative correction factor per se uncalculable in principle either 
with an accuracy ~ Aqcd/^6- In contrast, if one uses the short- distance masses, 
the higher-order corrections become smaller and the overall factor is theoretically 



calculable with the necessary precision [|l^, [22]. 



This, at first glance rather academic, observation in reality proved to underlie 
the pattern of the corrections from the very first terms. Remarkably, the actual 
model-independent calculations of F B i through observables measured in experiment 
are very stable against perturbative corrections. Including 0{a 2 s ) terms in the ex- 
traction of the b pole mass from, say, the e + e~ — > bb threshold region [2^] noticeably 



increases its value. However, if one performs the parallel perturbative improvement 
in calculating the width, one finds an essential suppression of the perturbative factor 
se pert . The two effects offset each other almost completely. 

This conspiracy is not unexpected: the appearance of large corrections at both 
stages is merely an artefact of using the ill-defined pole mass in the intermediate cal- 
culations. The situation is peculiar since the actual nonperturbative effects appear 
only at the level l/rriQ, whereas the pole mass is infrared ill-defined already at an 



accuracy of 1/niQ. The failure to realize this fact led to a superficial suggestion [24 
that even in beauty particles the perturbative corrections may go out of theoretical 
control; a more careful analysis |L5], showed that this is not the case. 

Since the apparent troubles with the perturbative corrections are associated with 
the pole masses, it is better to get rid of them altogether using instead running 
masses m(, )C (jtt) normalized at \i ~ 1 GeV; this is advantageous both theoretically 
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and in practice. This has been done in |T3[ and showed that neither masses nor 



the perturbative corrections to the width as pert have significant contributions from 
higher orders. It is important to note that it is not a kind of empirically selected 
procedure but is required in the literal implementation of Wilson's OPE in the case 
of the heavy quark widths [T9 |. 

To summarize, the idea that the perturbative corrections in the extraction of 
\V c b\ from T s i(B) are large emerges from an attempt to use ill-defined pole masses 
in an inconsistent way where two problems are encountered: 

• It is 'difficult to extract' accurately m^ ole from experiment; in any particular 
calculation the effects are easily identified that were left out, which can change 
its value by ~ 200 MeV. This uncertainty leads to a 'theoretical error' 5j in 
T A (B) of about 10%. 

• When routinely calculating r s i(i?) in terms of the pole masses, there are sig- 
nificant higher order corrections Sn ~ 10%. 



The naive conclusion drawn from such experience 
calculate the width without ~ 20% uncertainty: 



5r s 



6i + 5 U ~ 20% 



S\V C 



cb\ 



W, 



is that one cannot reliably 



10% 



cb\ 



On the contrary, theory predicts a cancellation between 5i and 5u in a consistent 
perturbative calculation, and that was explicitly checked in ]T^, [25]. For example, 
the net impact of the calculated (presumably dominant) second-order 0(a 2 s ) correc- 
tions on the value of \V c b\ appeared to be less than 1% ! Moreover, just neglecting all 
perturbative corrections altogether, both in the semileptonic width and in extracting 
nib from experiment, yields a \V c b\ smaller by less than 5% f27 |. 

Recently, all-order corrections associated with the effect of running of a s in one- 
loop diagrams (I shall refer to it in what follows as an extended BLM, or, shorter, 
merely a BLM approximation) were calculated in f25|, thus improving the exact 



one-loop result which has been known from the QED calculations [28 





28 


. Using 


23 


one gets 



most accurate model-independent low-energy determination of mb p3[ one gets fl5 

(BR{B 



cb\ 



0.0413 



0.105 



0.012 



(g - 0.4 GeV 
0.1 GeV 2 




(14) 



The relevant sources of theoretical uncertainty are shown explicitly. The main one is 
lack of knowledge of the exact value of /i 2 (marked with a tilde in eq. (|TJp, indicating 
that a particular field-theoretic definition of the kinetic operator has been assumed), 
which enters through the value of m^ — m c , eq. (0). A dedicated analysis of the lepton 
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spectra will reduce this uncertainty. At the moment a reasonable estimate of the 
actual uncertainty in is about 0.2 GeV 2 , leading to a 2.5% uncertainty in \y c \\. 

The dependence on the absolute value of m^ is minor; since we rely here on the 
well-defined short- distance mass ml, there is no intrinsic uncertainty in its value. 
The analysis |2B| estimated 5ml ~ 30MeV which seems reasonable; to be confident, 



I double it, and even then the related uncertainty in \V c b\ is less than 1.5%. 

Finally, we should consider the perturbative corrections. As previously ex- 
plained, the actual impact of the known perturbative corrections when relating the 
semileptonic width to other low-energy observables is very moderate, and there is no 
reason to expect the higher-order effects to be significant. With the a priori dom- 



inant all-order BLM corrections calculated [25[], one may be concerned only with 
the true two-loop effects 0(af) which are not associated with the running of a s . 



These have not been calculated completely yet; however, the recent calculation [2j| 
of the complete 0(a 2 s ) corrections in the small velocity kinematics demonstrated 
the feasibility of a complete calculation and suggested that the corrections must be 
small. There are some enhanced higher-order corrections not related to running of 
a s , that are specific to the inclusive widths PP[. They have been accounted for in 



the analyses (IS] , but went beyond those in |y, ^5J. With the results of |29 
it seems unlikely that as yet uncalculated second-order corrections can change the 
width by more than 2-3%; therefore, assigning an additional uncertainty of 2% in 
\V c b\ is a quite conservative estimate. 

Adding up these uncertainties we conclude that at present we can confidently 
calculate \V c b\ from T s \ with the theoretical uncertainty 



S\Vr 



cb\ 



< 5% . (15) 



c6 



Keeping in mind some controversy about experimental value of BR s i(5) and certain 
changes in time in tb, it is fair to say that the theoretical accuracy in extracting 
\Vcb\ is even better than its current experimental counterpart. 

It is important that this method has a potential for further improvement in 
a model-independent way. I think that the 2% level of a defensible theoretical 
precision can be ultimately reached here. It is difficult to count on an essential 
improvement beyond that, because of impact of higher-order power corrections and 
possible violations of duality. 

In a similar way it is straightforward to relate the value of \V u b\ to the total 
semileptonic width T(B — > X u tv) (TF 



llfcl- 0.00458 ( mB n Z; jW Wlgg V , (16) 



V 0.002 J V t b 

An accurate measurement of the inclusive b — > u£u width is difficult and for a 
long time seemed unfeasible. However, recently ALEPH announced the first direct 
measurement [j52 |: 

BR(B -> XJu) = 0.0016 ± 0.0004 . 
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I cannot judge the reliability of the quoted error bars in this complicated analysis; 
it certainly will be clarified in the future. Accepting the numbers literally, I arrive 
at the model-independent result 



IK 



ab | 



w, 



0.098 ±0.013 . (17) 



cb\ 



The theoretical uncertainty in translating T(B — > X u £u) into \V u b\ is a few times 
smaller and is not seen in the final result. 

To conclude this section on the inclusive semileptonic widths, let me briefly com- 
ment on the literature. It is sometimes stated that the existing uncertainty in T s \ 
is at least 20% [|26], |33|. The major origin of such claims is ignoring the subtleties 
related to using the pole mass in the calculations, and considering separately the 
perturbative corrections to the pole masses and to the widths expressed in terms 
of ^Q° le - This is inconsistent on theoretical grounds [pi|, whose relevance was con- 



firmed by the concrete numerical evaluations of the BLM corrections in [|15| , |25 
The dependence on nit, and nib — m c used to determine the uncertainty in \V C b\ was 
calculated erroneously in |2(| (cf. ]T5fl), apparently because of an arithmetic mistake 
that led to a significant overestimate. Finally, no convincing argument was given 
to justify a sevenfold boosting of the theoretical uncertainty in mj obtained in the 
dedicated analysis p3f . 



2.2 Exclusive zero-recoil B — > D* Iv rate 

Good theoretical control of all QCD effects in the inclusive widths was due to the fact 
that removing constraints on the final state to which decay partons can hadronize, 
makes such a probability a truly short- distance quantity amenable to a direct OPE 
expansion. Application of a similar idea to the exclusive zero-recoil decay rate B — > 
D* tv yielded quite an accurate determination of \V c b\ as well || [T7j], though with 
a more significant irreducible model dependence and a larger intrinsic uncertainty. 
The limitation is twofold: first, constraining the decays to a specific final state 
makes the transition not a genuinely short-distance effect; second, it suffers from 
downgrading the expansion parameter, namely l/m c rather than I /nib. 

On the quark level the decay is shown in Fig. la; near zero recoil it is determined 
by the single hadronic formfactor F^,*. In the absence of corrections violating the 
heavy quark symmetry, F D * = 1 holds; for finite m h c it acquires corrections: 

Fd* = 1 — (1 — Va) + 6 1/m ? + ... . (18) 

The effect of the nonperturbative domain denoted by the last term and the ellipses 
starts with the terms ~ l/(m c ,nib) 2 f34] , 35|, but otherwise is rather arbitrary, since 



it depends on the details of the long-distance hadronization dynamics in the form of 
wavefunction overlap. This opened the field for speculations and certain controversy 
about the value of F D * [[36| . 
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The situation as it existed by 1994 was summarized in the review lectures [57| : 

r] A = 0.986 ± 0.006 S 1/m 2 = (-2 ± 1)% , (19) 

yielding Fp* ~ 0.97, and was assigned the status of "one of the most important and, 
certainly, most precise predictions of HQET" . Nowadays we believe that the actual 
corrections to the symmetry limit are larger, and the central theoretical value lies 
rather closer to 0.9 g |T7||. 



Regarding the perturbative calculations per se, it was later pointed out f38 
that the improvement []39] of the original one- loop estimate was incorrect, and the 
proper value is rather tja ~ 0.965±0.025; subsequent calculations of the higher-order 
corrections in the BLM approximation |E5[ |40| confirmed it: rjA ~ 0.965 ±0.02. The 
purely perturbative chapter was closed recently with the complete two-loop 0(a 2 s ) 
result [f23| 

^(2ioo P ) = o.960 ±0.007 . (20) 
It should be noted, however, that the inherent irreducible uncertainty of the complete 



perturbative series for t]a exceeds the one stated in eq. ( |20"D by a factor of three 

11, im. 

If the mass of the charm quark were a few times larger, from the practical 
viewpoint the two-loop calculation would have been the whole story for Fp*. In 
reality, the power corrections originating from the domain of momenta below ~ 
0.6 GeV appear to be more significant. Unfortunately, not much can be said about 
them in a model-independent way, although they have been shown to be negative 
and exceed about 0.04 [§, [L7| in magnitude (the analysis for A b was given in |f44| ). 

The idea of this dynamical approach was to consider the sum over all possible 
hadronic states in the zero- recoil kinematics but not to limit it to only B — > D*; 
such a rate sets an upper bound for the production of D*. This inclusive quantity 
is of a short distance nature and can be calculated in QCD by means of the OPE. 
Schematically, the result through order 1/m 2 is 

\Fo*\ 2 +Em 2 = Uri-^ 2 -^^(- 2 +^ + ^) , (21) 

im A c 4 \vci c mi 6m c mb J 

where Fi are the-axial current transition formfactors to excited charm states i with 
the mass Mi = Mr>* ± e^, and £a is a perturbative renormalization factor (the role of 
will be addressed later). Considering a similar relation for another type of 'weak 
current', say, ci^b, one obtains a different sum rule 

H\h?={---) 2 ^^ , (22) 

with the tilde referring to the quantities for the transitions induced by this hypothet- 
ical current. These sum rules (and similar ones at arbitrary momentum transfer), 
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established in Refs. || [L7|], have been subjected to a critical scrutiny for two years, 
but are now accepted and constitute the basis for currently used estimates of F D * . 

Since eq. ( |2"2"D is the sum of certain transition probabilities, it is definite-positive 
and results in a rigorous lower bound 

/4 > /4~0.4GeV 2 • (23) 

The sum rule ( pT|) then leads to the model-independent lower bound for the 1/m 2 
terms in Fd*\ 

Ml* -M 2 

- 5 1/m2 > B ~ 0.035 . 24 

Smf. 

The actual estimate depends essentially on the value of /i 2 . It was suggested in || 
to estimate the contribution of the excited states in the l.h.s. of the sum rule (j2~ID 
from to 100% of the power corrections in the r.h.s.: 



J Ml* -Ml ul- u% (I 1 2 \\ 

_a 1/mJ = (1 + x) (^^ + U(_ + _ + _jj , < X <1. 

(25) 

If so, one arrives at |8j 

-5 1/rn 2 = (5.5 ± 1.8)% at /i 2 = 0.4 GeV 2 

- 5 1/m 2 = (6.8 ± 2.3)% at /i 2 = 0.5 GeV 2 (26) 

-8 1/m 2 = (8.1 ± 2.7)% at /i 2 = 0.6 GeV 2 

Before proceeding to further phenomenological implications, let me dwell on the 
QM meaning of the sum rules flT7| |. The act of a weak semileptonic decay of the 
b quark is its instantaneous replacement by c quark. In ordinary QM the overall 
probability of the produced state to hadronize to some final state is exactly unity, 
which (neglecting radiative corrections) is the first, leading term in the r.h.s. of (pi]). 
Why then are there nonperturbative corrections in the sum rule? The answer is that 
the normalization of the weak current 07^756 is not exactly unity and depends, in 
particular, on the external gluon field. Expressing the QCD current in terms of the 
nonrelativistic fields used in QM one has |T7|, E|, for example, 

1 1 _-, 2 1 ^.^ „/ 1 



C7fc7 5 & <-> -—r(aiD) a k + —^a k (aiD) - (aiD)a k (aiD) )+0 — 

\om^ om^ Am c mb ) \mr 

(27) 

The last term just yields the correction seen in the rhs of the sum rule. 

The inequality /x 2 > Hq in QM expresses the posit ivity of the Pauli Hamiltonian 

J_(^ l Df = — ((iD) 2 --aG 
2m v ' 2m V ' 2 

0. It is interpreted as the Landau precession of a charged (colored) particle in 
the (chro mo) magnetic field where one has (p 2 ) > \B\. Literally, in B* the QM 
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expectation value of the chromomagnetic field is suppressed, only /Xg/3, and in B 
it vanishes, since B is proportional to the spin of the light degrees of freedom. 
However, essentially non-classical nature of B (e.g., (B 2 ) > 3(B) 2 ) in turn enhances 
the bound which then takes the same form as in the external classical field. 

Before returning to numbers, let me emphasize that the perturbative factor £4 
is not (and cannot be) equal to r]\ |HJ (see also []19| ). It is clear that purely 
perturbative calculations include (although in an improper way) the effect of the 
strong interaction domain completely described by the power terms. £4 depends 
explicitly on the separation scale \x dividing the two domains. Similarly, the values 
of n 2 and \j? g depend, in fact, on \i (which, for simplicity, was neglected above). 

Most important, unlike r\ A which in principle cannot be defined theoretically 
with better than a few percent accuracy, £a(aO is a well-defined quantity and can be 
calculated in the small-coupling expansion provided \x is large enough. It is fair to 
say that no significant uncertainty can be associated with this factor which is close 
to unity, £ A ~ (0.99) 2 . 

What can be concluded for Fd* numerically? Even allowing a very moderate 
interval for the yet unknown value of /x 2 between 0.4 GeV 2 and 0.6 GeV 2 , we see 
—5i/ m 2 varying between 3.5% and 11%; moreover, since there are hardly any model- 
independent arguments to prefer any part of the interval a priori, the whole range 
must be considered equally possible. Adding small perturbative corrections we end 
up with the reasonable estimate Fjj* ~ 0.9 . It is curious to note that at a 'central' 
value x — 0.5 the dependence of the zero-recoil decay rate on /z 2 through 5y m 2 
practically coincides with that of T S \(B), see eq. (PJ), although they actually change 
in the opposite directions. 

The typical size of the l/m^ corrections to the exclusive zero-recoil decay rate 
seems to be significant, around 15%, which is expected since they are driven by the 
scale m c ~ 1.3 GeV. It is evident that 1/m;? correction in Fd* not addressed so far 
are at least about ±(0.15) 3 / 2 ~ 2-3%.[| 

Thus, I believe that the current theoretical technologies do not allow to reliably 
predict the zero-recoil formfactor Fd* with a precision better than 5-7% in a model- 
independent way; its value is expected to be approximately 0.9, although a correction 
to the symmetry limit twice smaller, as well as larger deviations, are possible. It 
is encouraging that the 'educated guess' Fd* ~ 0.9 which emerged from the first - 
and so far the only - dynamical QCD-based consideration JSL [Tj]], yielded a value 
of \V c b\ rather close to a less uncertain result obtained from T s i(B), although the 
experimental error bars for the zero-recoil rate are still significant and its status 
does not seem finally settled yet. 

In future, more accurate data will enable us to measure Fd* with a theoreti- 
cally informative precision using \V c b\ from T S \(B), and thus provide us with deeper 
insights into the dynamics of strong forces in the heavy quark system. 

1 This is consistent with the fact that the 1/tUq IR renormalon ambiguity in r)\ constitutes 5% 
at A^ D ~ 220 MeV g§. 
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Smaller theoretical uncertainty, 2.5% and 3%, is now quoted by Neubert for 5jy m 2 
and Fd*, respectively. The former was obtained in ||4lf , in what he calls a "hybrid 
approach" ,[] which reduces to assigning the fixed value fj, 2 = 0.4 GeV 2 and using it 
in the sum rule (21) within the same model assumption < x < 1 suggested in 
M. Correspondingly, the quoted number for 5i/ m 2 practically coincided with the 
first line of eqs. (^). In reality, allowing fi 2 to vary within any reasonable interval 
significantly stretches the uncertainty. Moreover, the analysis [[|l|] was based on 
using t]a as the perturbative factor assuming, literally, that in the proper treatment 
the final result would not be changed numerically - which is just the case according 
to consideration of Ref. 
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On top of that, the uncertainty in the definition of 
t]a due to 1/m 2 and 1/m 3 IR renormalons constitutes 2-3% each and thus is an 
additional one in the usage of the sum rules adopted in pEIfl. As a result, the stated 



theoretical accuracy of those estimates cannot be accepted as realistic. 



3 Perturbative corrections to the sum rules and 
the IR renormalons 

The perturbative corrections have been mentioned already in the context of the 
semileptonic widths; the nonperturbative effects there were small. A different situa- 
tion occurs when the latter are in the focus; the interplay between the perturbative 
and nonperturbative physics becomes important. This happens in the analysis of 
the sum rules, in particular, when the effects of the running of a s is included. 

A typical short- distance observable A at the one-loop level can be schematically 
written in the form 

MQ') ~ i + /0^)«-(* 2 > (») 

where <fi is some function and Q 2 sets the high momentum scale. Even at Q 2 ^> 
Aq CD , small momenta k 2 ~ Aq CD ~ Q 2 e _47r /( 6a s) contribute at a certain level, but 
in this nonperturbative domain the coupling a s (k 2 ) cannot be expanded in terms 
of a s (Q 2 ). This leads to emergence of the factorial growth of the coefficients in the 
perturbative expansion of A: 

A{Q 2 ) ~ 1 + f] a n , a n ~ ^ n\ . (29) 

This (same-sign) factorial behavior manifests the presence of the IR renormalon 
singularity |4^] which makes the series ill-defined at the level of power corrections, 

2 I disagree with the statements of EJ, reiterated in later papers, suggesting that the original 
analysis ^ [Tt]] missed some elements of the heavy quark spin-flavor symmetry; on the contrary, 
it was stated in the latter paper that all these relations automatically emerge from the sum rules 
that replace the QM wavefunction description in the quantum field theory. 
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~ (Aqcd/Q)', I = 1, 2, ... . On the other hand, this is a deficiency of the purely per- 
turbative expansion itself and has nothing to do with the actual strong-interaction 
domain [|47||; in particular, the factorial growth persists even if the effective cou- 



pling never becomes large but stays finite for all k 2 , when the integral in eq. (p8| ) is 
unambiguously defined E7L E8 fl. 



It is important to note that the IR renormalons are absent in the Wilson OPE, 
when the separation between the operators and the coefficient functions is done ac- 
cording to the short-distance or long-distance origin of the contribution |J. They 
appear only at the attempts to define 'purely perturbative' and 'purely nonpertur- 
bative' pieces of an observable. In this respect processes with the heavy quarks do 
not have any peculiarity. 

The IR renormalons in the HQET were addressed recently also in Refs. E3, HU 



I cannot agree with the message of ||9 which claimed the IR renormalons to be 



a peculiar feature of the effective theories. On the contrary, the IR renormalons 
are generally identical in QCD and its effective low-energy limit, since their differ- 
ence, by definition, lies only in the ultraviolet domain. The structure of the IR 
singularities in perturbative QCD does not depend on whether the QCD is a fun- 
damental underlying theory, or is an effective low-energy representation of a more 
general quantum system like GUT or a superstring theory. Appearance of the IR 
renormalons in the perturbative calculations depends solely on the way how one 
performs the OPE in a theory. 

A certain peculiarity of the HQET in its standard formulation, in this respect, 
is that it has an extra spurious 1/toq IR renormalon. It does not have any deep 
reason but is merely rooted in the definition of the field variable routinely used in 
the HQET for the heavy quarks (for simplicity I ignore the spinor indices): 

h v (x) = e i{vx)m Q 1C -Q(x) (30) 

with the pole mass in the exponent. It is just this factor that brings in the leading 
1 / tuq renormalon singularity into the HQET perturbative series for short- distance 
quantities. Since none of the relations between the observables can depend on the 
phase definition used in the calculations, this spurious effect always disappears in 
the relations among short-distance observables At the technical level, however, 
it sometimes may look not obvious. 

In the Wilson OPE one defines an effective theory normalized at some scale //; in 
particular, the basic parameter for the heavy quark physics is mg(/i), which is the 
short-distance quantity, not sensitive to the physics below fi. The most convenient 
(but, certainly, not the only possible) choice for the nonrelativistic field definition is 
then 

Q M (x) = e*^ m ^-Q M (x) (31) 
where the subscript (/x) denotes the normalization point. In this case one has an 
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operator equation of motion 




t(vD) M Q M = + O . (32) 



In the perturbation theory the product of the operators in the l.h.s. has the pertur- 
bative corrections ~ a" • \x coming from both the operator D a = d a + g s A a and from 
the interaction of Q with the gluon field. They combine to satisfy eq. (|32|) to any 



order |19| . In the standard HQET approach one insists (4^, [EJ that no perturba- 
tive correction is present from any separate source, and this unnatural requirement 
leads to the IR renormalons in the matching with the QCD. 

Based on the renormalon analysis, Neubert claimed [fl], fEJ that the sum rule 
(PT|) cannot be correct, since I/w/q renormalons allegedly mismatch in it. Although 
in Wilson's OPE the IR renormalons are always absent from any particular term, 
the IR renormalon calculus can still be applied if the OPE relation is considered in 
the pure perturbation theory itself, and formally setting /i = (which is technically 
possible to any finite order in a s ). This is nothing but a cross-check how the OPE 
works in the one-loop perturbative calculations. Compared to Wilson's procedure, 
the limit // — ► 0, in particular, amounts to subtracting a 'perturbative piece' from 
the observable probabilities entering, for example, the phenomenological side of the 
sum rules. In this way the perturbative terms obviously appear in the left-hand 



side of eq. ( |2"T| ) as well, and these terms were ignored in [Efl . Including them 



immediately balances - not surprisingly - the IR renormalons in the both sides [|3| 



On the other hand, the analysis of Ref. ||2| is instructive demonstrating once 



more that HQET in that interpretation is not a closed theory once nonperturbative 
effects are addressed. For example, the sum rules in their literal form indeed cannot 
be derived in it. Clearly, there is nothing wrong with the sum rules themselves, 
the reason is merely non-existence (at the nonperturbative level) of HQET as an 
effective quantum theory with //-independent matching coefficients, according to 
prescriptions in its popular formulation |36[ . 

Leaving aside theoretical aspects of the IR renormalons, it is still important to 
calculate the perturbative corrections in the practical applications of the OPE. The 
first such calculation for the sum rules was done in |T7[ and then addresses again 



in a number of papers fig] . The 0(a s ) terms appeared to be rather moderate. 
Based on the first terms ~ ha 2 s in the BLM series, it was suggested in [[|!| that 
the higher-order radiative corrections to the sum rules are too large and allegedly 
make them next to useless. These conclusions stemmed from the alternative, non- 
Wilsonian use of the sum rules of H, [I7j when the perturbative pieces as they literally 
emerge in the BLM calculation, are completely subtracted from all the entries. The 
large corrections then came basically from integrating the running coupling in the 
domain below the Landau singularity. The numerical analysis in the OPE |43|, on 
the contrary, suggested a quite moderate impact of the radiative corrections. 



According to p5l, the perturbative corrections to the sum rule of the type of 



eq. fl22|) weaken the bound for the expectation value of the kinetic operator to such 
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extent that it becomes non-informative. To appreciate the status of this statement, 
one must realize that, in the quantum field theory, the renormalized operators can 
be defined in different non-equivalent ways; — Ai addressed in is known to be 
different from /^ 2 . In general, the complete definition of an operator at the nonper- 
turbative level is a nontrivial task. As a matter of fact, the only complete definition 
of the kinetic operator Q(iD ) 2 Q given so far was made in JT^j. In any case, an effec- 
tive operator, intrinsically normalization- scale independent as intended in HQET, 
cannot be defined in a QCD-like theory. The inequality > \j? g always holds 



for the expectation value of the well-defined operator of [17|, for arbitrary normal- 



ization point. Moreover, with this definition the perturbative normalization-point 
dependence has been calculated there (see also [[43|1 ). 

As for — Ai, a parameter in HQET, its definition beyond the classical level has 



never been given; the procedure adopted in |45| reduces to an attempt to completely 
subtract the 'perturbative piece' of fJ%{fJ>)' 



a s (ji) 2 

Ci—fi 
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C'2 
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(33) 



(the method to calculate q was elaborated in fll7fl ). Yet such a program cannot be 
performed even in principle ||. The situation appears rather simple in the BLM 
approximation, where all q are readily calculated: the series, whose second term 
was discussed in |E| , is divergent and sign- varying , so using merely the second 
term is misleading for evaluation. Numerically, the perturbative subtraction looks 
as follows: 



-Ai =/4(/i) -0.12 GeV 2 -0.3GeV 2 + ... > 0.4 GeV 2 -0.12 GeV 2 -0.3 GeV 2 + ... , 

(34) 

with the ellipses denoting higher-order terms which already grow. Based on the 



third term, it was concluded in ]45] that no bound for the expectation value of the 
kinetic operator can be derived. In what concerns the HQET parameter — Ai, it is 
perfectly true: clearly there can be no definite bound established for a quantity that 
has not been - and cannot be - defined. 

Concluding the discussion of the kinetic operator, I should emphasize that the 
above subtleties are peculiar to the field-theory analysis. The inequality /i 2 > \Xq 
must hold in any QM model relying on a potential description not invoking ad- 
ditional degrees of freedom beyond the heavy quark and a spectator, if the heavy 
quark Hamiltonian is consistent with QCD. Unfortunately, a failure to realize this 
general fact is seen in a number of recent analyses. 

A similar second-order BLM improvement of the sum rule (pTf) for Fd* was also 



attempted in [45] and claimed to destroy the predictive power of this relation (the 
first-order calculation had been performed in |T7[ as well). The calculation of the 
impact of the second-order perturbative terms, however, was not done consistently, 
and the actual effect is smaller [43| . 
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Figure 3: The value of r] A (fi) for Aqq D = 90 MeV, 130 MeV and 175 MeV. The value 
of rjA {n) must be used in the QCD-based calculations of the exclusive zero-recoil 
B — > D* form factor when nonperturbative effects are addressed. The shaded bars 
show the perturbative uncertainty irreducible without considering \jm\ and l/m^ 
effects; they give the lower bound for the corresponding actual nonperturbative 
corrections. 



In reality, the perturbative corrections seem to be modest, at least for reasonable 
values of the running coupling. To illustrate this assertion without submerging into 
details, let me introduce the following notation: 

VaM = IU^)} 1/2 • (35) 

77a(aO can indeed be calculated in the perturbative expansion; for example, to order 
a s it is given by 



Va(p) = 1+ — 

7T 



rrib + m c nib 8 1 / /i 2 /i 2 2/i 



nib — m c m c 3 3 \ m 1 mi 3m c mb 



(36) 



The quantity t]a{ij) must be added to Sy m 2 in the framework of the OPE to calculate 
or bound Fd* , instead of t]a in the model calculations. Then, at a reasonable choice 
jj, ~ 0.5 GeV, Aqcd = 300 MeV (in the V-scheme) one has 

t}a(/J') = 1 tree level 

t)a{ia) = 0.975 one loop 
T] A (ii) = 0.99 all-order BLM 

Clearly, the effect of the calculated perturbative corrections is not significant and 
T]a{^) is very close to the value of 0.98 adopted in the original analysis 0. The 
dependence of the BLM-resummed value of r) A (/j,) on /x and on the QCD coupling is 
illustrated in Fig. 3 (Aqcd is shown in the MS scheme). 
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Conceptually, the deficiency of the alternative usage of the sum rules of Refs. || 
[L7j chosen in [[45| is to gauge £4 on the value of r] A as it has been defined in the 
HQET (the idea that £4 is to be identified with r)\ ascends to |^2], |^I|): the value of 



rjA is defined as an all-order 'purely perturbative' renormalization factor of the zero- 
recoil axial current. However, it is this factor tja that is ill-defined, and only for this 
reason must the difference between the stable Wilson coefficient £4 and n A suffer from 
large corrections. It is worth reiterating once more in this respect that, in reality, 
t\a cannot be equal to a matching coefficient of 07^756 to a corresponding current 
in any effective field theory completely defined at the level of the nonperturbative 
physics. 

The accurate analysis thus shows that the literal implementation of Wilson's 
OPE not only eliminates purely theoretical problems with the IR renormalons, 
but also makes the actual impact of the higher-order perturbative corrections well- 
controlled. 



4 Bounds on the transition formfactors from an- 
alyticity and unitarity 

The kinematic point of zero recoil in the b — > c semileptonic transitions is the least 
uncertain one for the theoretical analysis. On the other hand, experimentally the 
statistics is limited here. Therefore, information about the g 2 -dependence of the 
formfactors is very important. However, is has a dynamical origin and is unknown 
a priori. Indeed, if the B meson is relatively 'soft', i.e. the typical velocity of the 
effective light degrees of freedom is small - as would be the case in the nonrelativistic 
quark models - then the formfactor (the IW function) changes at the velocity transfer 
much less than unity, and this implies large derivatives dF(v 2 )/dv 2 , d 2 F(v 2 ) / (dv 2 ) 2 , 
etc. If, on the contrary, the bound state is relatively simple and the light degrees 
of freedom are ultrarelativistic, these derivatives are expected to be of order unity. 
Can something be stated about them without actual information about the bound 
state dynamics? I will argue below that the answer is negative. 

An idea was put forward a few years ago to constraint the IW function using, 
mainly, general analyticity and unitarity properties of the formfactors in the heavy 
quark limit They turned out to be non- informative due to the presence of the 



QQ-bound states below the threshold of the open flavor production. This approach 
was resurrected recently in a number of publications |52|, |K| . In particular, it was 
claimed in [53| to have established an intriguing relation between the slope and 



the curvature of the formfactor near the zero recoil kinematics. However, there are 
essential flaws in that analysis. 

The general idea is to consider the transition formfactor for some cb current 
between B and D mesons as an analytic function of the momentum transfer q 2 : 



(D(p + q)\ c6(0) \B(p)) = 2dM B M D F s (q 2 ) (37) 
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where, for simplicity, we can consider the scalar current not yielding irrelevant 
Lorentz indices. 

The function F$ in the physical t-channel domain q 2 < t zr = (Mb — Mo) 2 
describes the decay amplitude and is real. On the other hand, analytically continued 
to large q 2 > (Mb + Mo) 2 it describes the exclusive production of BD in e + e~- 
annihilation-type process and thus, in principle, is also an observable quantity. F$ 
is complex there. Under the usual assumptions of analyticity one can write the 
dispersion integral for F5: 



I r Im F s (t + ie) dt = 1 r lm F s (t) 

tt J t — q 2 — ie 7r Jt<(M B +M D ) 2 t — q 2 — ie 



+ - I ^ ; rL <^ = F s\t) +n'\t) ■ (38) 



7r Jt>(M B +M D ) 2 t — q 2 — ie 
If the sub-threshold contribution from lm.Fs(t) for t < (Mb + Mjj) 2 were small, 
the formfactor would be saturated by the second integral, F$ — Fg. In that case 
the formfactor would be essentially constrained: since the exclusive production of 
BD must not yield the production cross section in this channel larger than the total 
hadronic one which is bound theoretically, F cannot be too large in the wide domain. 
If one additionally has the normalization 

F s (t ZI ) = 1 (39) 

(which holds at large ra& )C ) then only little room would have been left and, for 
example, an almost functional relation between F s (t ZT ) and Fg(t zr ) emerged p3 . 



This reasoning, however, crucially relies on the insignificance of the sub-threshold 
support of ImF(g 2 ); on the other hand, it is generally known not to hold for heavy 
quarks. On the contrary, the common wisdom is that the first resonances dominate 
the formfactor. The idea of Ref. was to consider the formfactor Fg which does 



not have a contribution from the well-studied vector quarkonia. It was assumed that 
no clear-cut resonances below the BD threshold exist in the corresponding channel. 

This basic assumption is erroneous, however. The sub-threshold scalar ++ states 
are present both in the cc and bb systems (x c ,b states) and, according to the general 
theorems, must be also in the be channel.^ They are extremely narrow since lie 
below the open-flavor thresholds, and are expected to dominate the formfactor.^ 

To account for the possible effect of the domain q 2 < (Mb + Mo) 2 , Ref. |53 



equated the whole sub-threshold contribution with the possible nonresonant one 
and considered a smooth model for Im F$ below the threshold: 



lmF s (t) ~ f|== ^= ^t-(M Bc + ^) 2 , (M Bc +m v ) 2 < t < (M B +M D ) 2 

(40) 



3 In any case it is ensured by the heavy flavor symmetry which is inherent in the analysis of |53| . 

4 It can be argued that the contribution of the continuum domain q 2 > (Mb + Mo) 2 itself falls 

.2 _ c\/t m\2 



short of 1 for F$ at the zero recoil point q- = (Mb — Mb) 
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with C of order 1. As explained above, such a functional form is irrelevant below 
the threshold. Moreover, it is easy to calculate its contribution to the zero-recoil 
value F S ((M B -M D ) 2 ) : 

i^ r) (0) ~ 0.008 C (41) 

(recall that in the usual pole- dominance models this value must be just close to 
unity). Thus, the effect of the states lying below the open flavor threshold was 
underestimated by two orders of magnitude. 

The analysis of Ref. |k| applied to the scalar cb formfactor. Obtaining model- 
independent information about any formfactor would be very interesting by itself 
from the theoretical viewpoint. Phenomenologically, however, one needs to know the 
axial- vector formfactors, or those in the heavy quark limit itlq — > oo. In the heavy 
quark limit all these formfactors are related to the IW function; unfortunately, the 
l/m c corrections in the derivatives are too significant. Although they were claimed 
to be very mild, in reality it was based only on some perturbative estimates. On 
the other hand, the main, nonperturbative l/m c corrections are known to be large 

g,~5o%.g 

Can these problems be cured? One could have tried to consider similar dispersion 
relations at arbitrary masses m Ci &. Decreasing m c would soften the problem of the 
sub-threshold physics - however, the corrections in the short-distance expansions and 
to the necessary heavy quark symmetry relations are already marginal. Attempts 
to increase m c to tackle the l/m c corrections magnify the role of the sub-threshold 
behavior of the formfactor. 

It is important that, contrary to the forward scattering amplitudes or current 
correlators, the phase of the formfactor has only an indirect physical meaning, and 
below the threshold is a very unphysical quantity not observable directly. For this 
reason its behavior in this domain can be rather odd and crucially depends on 
irrelevant details of the strong interactions which we know almost nothing about; 
by no means, for example, Im F must be positive for heavy quarks. 

A very instructive analysis had been given a few years ago in the inspiring papers 
|H| which, unfortunately, are to a large extent forgotten in the recent literature. It 
was argued there that, just if the heavy flavor transition formfactor is governed by 
relatively soft clouds of light degrees of freedom, in the s-channel one should expect a 
complicated oscillating behavior of the production formfactor near the threshold; it 
can be simple if the transition formfactor is a slowly varying function (a 'hard-core' 
bound state) and then the lowest-pole dominance is expected to work. 

One has to conclude, once again, that without actual dynamical input, the gen- 
eral analyticity and unitarity relations cannot per se yield useful information about 
the transition formfactors of interest. From the phenomenological perspective, the 
analysis |k| grossly underestimated both the sub-threshold contributions and the 
corrections [^4] to the heavy quark symmetry constraints it relied on. The relation 

5 Although it was admitted in that the impact of the l/m c corrections was missed in the 
analysis, this qualification was somehow absent from the phenomenological claims. 
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between the slope and the curvature claimed there rather should not be used in 
experimental analyses for deriving model-independent results. 

The possibility to employ additional dynamical information, e.g. in the form of 
the lattice calculations of the formfactors at some intermediate points [p6j, must be 
explored in order to revive this approach. 



5 The semileptonic branching ratio 

In the QCD-based heavy quark expansion the total nonleptonic decay widths of 
heavy flavors can be treated in the same way as semileptonic or radiative ones; the 
difference appears only at a quantitative level. The statement about the absence 
of 1/rriQ corrections applies as well, and the leading corrections ~ 1/^q have been 
calculated || [7|, |12| . The overall semileptonic branching ratio BR s i(5) seems to be 
of a particular practical interest: while the simple-minded parton estimates yield 
BR sl (£?) ~ 15% |57| , experiment gives smaller values, BR sl (i?) ~ 10.5-11.5% (the 



current situation is reviewed by T. Browder). 

The generic scale of the 1 / m\ effects in the widths of individual nonleptonic chan- 
nels is about 15%; however, the particular V-A chiral and color structure of weak 
currents significantly suppresses it; moreover, an additional cancellation between 
different parton channels occurs, and one literally gets |58| 



A 1/m2 (BR sl ) £ 0.5% . (42) 



Estimated \jm\ corrections do not produce a significant effect either [59]. As a 
result, most of the attention was paid to a more accurate treatment of the pertur- 
bative corrections, in particular, accounting for the effect of the charm mass in the 
final state. It was found |ET| that the nonleptonic width is indeed boosted up. 

Since the inclusive widths are expanded in inverse powers of energy release, one a 
priori expects larger corrections or even a breakdown of the expansion and violation 
of duality in the channel b — > ccs(d); however, this channel can be isolated |J via 
charm counting, i.e. measuring the average number of charm quarks per beauty 
decay. The original experimental estimate n c ^ 1.15 did not allow one to attribute 
the apparent discrepancy to this channel, and gave rise to the so-called 'BR sl versus 
ric problem. 

The perturbative corrections in the b — > cud itself cannot naturally drive BR s i 
below 12.5%; the calculation for the b — > ccs(d) channel is less certain and, in 
principle, admits increasing the width by a factor of 1.5 or even 2, leading to n c ~ 
1.25-1.3. In the latter case a value of BR s i as low as 11.5% can be accommodated. 

The experimental situation with n c does not seem to be quite settled yet: 

n c = 1.134 ±0.043 (CLEO) , n c = 1.23 ±0.07 (ALEPH) . 



In a recent analysis ]62| it was argued that consistency requires a major portion 
of the final states in the b —>■ ccs to appear as modes with the wrong-sign .D's 
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and kaons but D s , which previously escaped proper attention, and this allowed for 
a larger value n c ~ 1.3 needed to resolve the problem with BR sl . The dedicated 



theoretical consideration shows that, indeed, the dominance of such modes is 
natural and this possibility does not require a significant violation of duality. Thus, 
if a larger value of n c 1.25 is confirmed experimentally, the problem of BR s i will 
not remain. 

Even within this successful scenario a concern can be raised if a strong QCD- 
enhancement of a tree-level-unsuppressed decay mode, possibly achieved in the one- 
loop calculations, is trustworthy: one would need to boost the tree-level b — > ccs 
width by a factor of two. One should note, however, that theoretically the situation 



maybe not so pessimistic: the tree level estimate of BR s i(6 — > ccs) ~ 15% [57 
refers only to using the large values of the quark masses close to their pole values. 
Using the short-distance masses like m& — 4.5 GeV always yielded a larger starting 
value of the branching for this channel. Without addressing the strong interaction 
effects it must be viewed, conservatively, as an uncertainty of the simplest partonic 
predictions. However, it is more consistent, from the viewpoint of the OPE, to always 
use the values of the parameters entering at a particular level of computations; in the 
case of the B widths it means that the smaller short-distance masses are preferable 
as the starting approximation. If one accepts this attitude, the actual corrections 
to the decay rate found at the one-loop level are rather mild and do not prompt 
immediate concerns about the theoretical control over the perturbative expansion. 
The observed large corrections to BR s i(6 — > ccs) are then merely a response to an 
improper zeroth-order approximation. 



6 Lifetimes of beauty particles 

The Q CD-based HQE provides a systematic framework for calculating the total 
widths of heavy flavors, which are not amenable to the traditional methods of HQET. 
The only assumption is that the mass of decaying quark (actually, the energy release) 
is sufficiently large; for a review, see |64 . 



A thoughtful application of this expansion to charm particles demonstrated that 
the actual expansion parameter appeared to be too low to ensure a trustworthy 
accurate description, so that a priori one expects only emergence of the qualitative 
features. Surprisingly, in many cases the expansion works well enough even numer- 
ically. In particular, it predicts « 2t£,o as an effect of Pauli interference in 
D + and suggests Tn a t d o - simultaneously correctly predicting the semileptonic 
fraction in D mesons. A similar agreement seems to be observed in charm baryons 
(for a recent review, see |JS| ) . 



Applying the expansion to beauty particles one expects a decent numerical accu- 
racy, although the overall scale of the preasymptotic effects is predicted to be small 



24 



and makes a challenge to experiment: 

t b -/t b o ~ l + 0.04( T g 5 % v ) 2 [59] EXP: 1.04 ±0.04 

t b Jt b o ~ 1 + 0(1%) @ EXP: 0.97 ±0.05 



{T Bi -TBs)/T Ba * 0.18 (55^) II 



(43) 



T-A ft /r B o « 0.9 EXP: 0.78 ±0.06 



These differences appear mainly as l/rn| corrections and, depending on certain 
four-fermion matrix elements, cannot be predicted at present very accurately. In 
particular, it refers to the preasymptotic effects in baryons. For mesons the estimates 
are based on the vacuum saturation approximation, which cannot be exact either. 
The impact of non-factorizable terms has been studied a few years ago in ]67|] and 



possibilities to directly measure the matrix elements in future experiments were 
suggested. 

The apparent agreement with experiment is obscured by the reported lower val- 
ues of ta 6 . Since the baryonic matrix elements are rather uncertain, quite a few 



model estimates have been done (see the recent one |]68 |, and references therein). 
All seem to fall short; however, this might be attributed to deficiencies of the simple 
quark model. It was shown using quite general arguments that, irrespective 
of the details, one cannot have an effect exceeding 10 — 12% while residing in the 
domain of validity of the standard I/tuq expansion itself; the natural 'maximal' 
effects that can be accommodated are ~ 7% and ~ 3% for weak scattering and 
interference, respectively. 

Thus, if the low experimental value of ta 6 is confirmed, it will require a certain 
revision of the standard picture of the heavy hadrons and of convergence of the 
l/m,Q expansion for nonleptonic widths. 



7 expansion and duality violation 

The problem of duality violation attracts more and more attention of those who 



study the heavy quark theory; a recent extensive discussion was given in | 70| . The 
expansion in I/uiq is asymptotic. There are basically two questions one can ask 
here: what is the onset of duality, i.e. when does the expansion start to work? 



The most straightforward approach was undertaken in Ref. JTTJ] , and no apparent 
indication toward an increased energy scale was found. Another question, how is the 
equality of the QCD parton-based predictions with the actual decay rates achieved, 
was barely addressed. Though a relevant example of such a problem is easy to give. 

The OPE in full QCD ensures that no terms ~ l/m^ can be in the widths and 
the corrections start with l/nig. However, the OPE per se cannot forbid a scenario 
where, for instance, 

6T h q r sin (m Q p) x 
^ G , p ~ A QCD . (44) 
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In the actual strong interaction, m& and m c are fixed and not free parameters, so, 
from the practical viewpoint these types of corrections are not too different - but the 
difference is profound in the theoretical description! It reflects specifics of the OPE 
in Minkowski space, and such effects can hardly be addressed, for example, in the 
lattice simulations. Their control requires a deeper understanding of the underlying 
QCD dynamics beyond the knowledge of first few nonperturbative condensates. 

In fact, the literal corrections of the type of eq. (|4]) are hardly possible; the 
power of 1/ rriQ in realistic scenarios is larger, and these duality corrections must be 



eventually exponentially suppressed though, probably, starting at a higher scale [70 



But a theory of such effects is still in its embryonic stage and needs an additional 
experimental input as well. 



8 Conclusions 

I would like to conclude this admittedly incomplete review of the theory and ap- 
plications of the heavy quark expansion by the following general remark. It was 
sometimes said already a few years ago that the theory of heavy quark decays was 
basically completed and only a few technical refinements remained to be done. The 
subsequent development clearly showed that the potential for further studies was 
vast. Regarding the current status, I believe that, although the basic principles of 
the HQE in QCD are formulated and well studied, we still have many interesting 
things to understand and calculate even in the original, the most standard fields of 
applications - before one can actually consider the theory completed. 
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